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We introduce a novel class of parametric optical solitons supported simultaneously by two second- 
order nonlinear cascading processes, second-harmonic generation and sum-frequency mixing. We 
obtain, analytically and numerically, the solutions for three-wave spatial solitons and show that the 
presence of an additional cascading mechanism can change dramatically the properties and stability 
of two-wave quadratic solitary waves. 



_j . As is known, optical cascaded nonlinearities due to parametric wave mixing can lead to a large nonlinear phase 
' shift and spatial solitary waves, resembling those for a Kerr medium However, solitary waves supported by 
^\ , cascaded nonlinearities demonstrate much richer dynamics due to nonintegrability of governing nonlinear equations 
C^I ' and, unlike solitons of the Kerr nonlinearity, the quadratic solitons can become unstable in a certain narrow region 
of their parameters [|j . 

" ^ ■ In this Letter we introduce a novel class of parametric spatial solitons supported simultaneously by two nonlinear 
QQ I quadratic (or x^^^) optical processes: second-harmonic generation (SHG) and sum- frequency mixing (SFM). As has 
been recently shown by Koynov and Saltiel for continuous waves, under the condition that the two wave-mixing 
. processes are nearly phase matched, the presence of multistep cascading leads to a four fold reduction of the input 

■ intensity required to achieve a large nonlinear phase shift. Here, we demonstrate that the multistep cascading can 
, lead to a new type of parametric solitons. Introducing a third wave generated via a SFM process, we find that it can 
' alter both the general properties and stability of the two-wave x*-^-* spatial solitons. Moreover, we reveal the existence 

. , of a new type of the so-called quasi-soliton, that appear for a negative mismatch of the SFM process. 
Q To introduce the model of multistep cascading, we consider the fundamental beam with frequency uj entering a 

, noncentrosymmetric nonlinear medium with a x^^-* response. As a first step, the second-harmonic wave with frequency 
■{— { ' 2a; is generated via the SHG process. As a second step, we expect the generation of higher order harmonics due to 
d , SFM, for example, a third harmonic {oj + 2uj ~ 3oj) or even fourth harmonic {2uj + 2lo — Auj) [Q. When both 
such processes are nearly phase matched, they can lead, via down-convertion, to a large nonlinear phase shift of the 
^ , fundamental wave Q. Additionally, as we demonstrate in this paper, the multistep cascading can support a novel 
type of three-wave spatial solitary waves in a diffractive x'^' nonlinear medium, multistep cascading solitons. 
We start our analysis with the reduced amplitude equations derived in the slowly varying envelope approximation 

■ with the assumption of zero absorption of all interacting waves (see, e.g., Ref. ||]). Introducing the effect of diffraction 
in a slab waveguide geometry, we obtain 

f) A ffi A 

2^fcl^ + ^ + Xi^3v4^e-^'=3- H- X2A,Ale-'^^'^-^-^ = 0, 

r) A rfl A 

4zfci^ + ^ + X4A3AJe-^'=3- + X5A?e^^'=^^ = 0, (1) 
6zfci^-H^ + X3A2Aie'^'=- = 0, 



where xi,2 — 2kiai^2, X3 — Gfeids, and X4,5 = 4fci(T4_5, and the nonlinear coupling coefficients ak are proportional to 
the elements of the second-order susceptibility tensor which we assume to satisiy the following relations (no dispersion), 
CTs = 3(71, a2 = cr^, and 174 = 2ai. 

In Eqs. (U), Ai,A2 and A3 are the complex electric field envelopes of the fundamental harmonic (FH), second 
harmonic (SH), and third harmonic (TH), respectively, Afc2 = 2ki — k2 is the wavevector mismatch for the SHG 
process, and Ak^ = fci -I- /c2 — fca is the wavevector mismatch for the SFM process. The subscripts '1' denote the FH 
wave, the subscripts '2' denote the SH wave, and the subscripts '3', the TH wave. Following the technique earlier 
employed in Refs. ||^, we look for stationary solutions of Eq. (|l|) and introduce the normalised envelope w{z,x), 
v{z, x), and u{z, x) according to the relations. 
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where A/c = AA;2 + Afcs. Renormalising the variables as z z/ (3 and a: ^ x/ \/2j3ki, we finally obtain a system of 
coupled equations, 

.dw d'^w ^ ^ 

i-^ — ^ -TT^ — w + w V + v u — I), 
oz ox'^ 

3«Tr + ~ "1" + Xvw = 0, 

where a = 2(2(3 + Ak2)/(3 and ai — 3(3/3 + A/c)//3 are two dimensionless parameters that characterise the nonlinear 
phase matching between the parametrically interacting waves. Dimensionless material parameter x = X1X3/X2 — 
9((Ti/ct2)^ depends on the type of phase matching, and it can take different values of order of one. For example, when 
both SHG and SFM are due to quasi-phase matching (QPM), we have aj — (2/7rm)(7r/Aini)x*^^) [w; (4— — (3— 
where j — 1,2. Then, for the first-order (m = 1) QPM processes (see, e.g., Ref. Q), we have cti = (T2, and therefore 
X = 9. When SFM is due to the third-order QPM process (see, e.g., Ref. 0), we should take ci = and therefore 

X = 1. At last, when SFM is the fifth-order QPM process, we have ai = 172/5 and x = 9/25. 

Dimensionless equations (|^) present a fundamental model for three- wave multistep cascading solitons in the absence 
of walk-off. Additionally to the type I SHG solitons (see, e.g., Refs [^), the multistep cascading solitons involve the 
phase-matched SFM interaction (lu + 2lu — 3uj) that generates a third harmonic wave. If this latter process is not 
phase-matched, we should consider ai as a large parameter, and then look for solutions of Eq. (|^) in the form of an 
asymptotic series in ai. Substituting w = wq + ewi + . . ., v = vq + svi + . . . and u — eui, where e = a^^, we find 
ui « x''Jw, and the system (||) reduces to a model of competing nonlinearities, 

.dw d'^w , I ,2 r> 

t-^ + — w + w V + exM w = 0, 

In the limit e — ^ 0, Eqs. (P coincide with the model of two-wave solitons due to the type I SHG earlier analysed in 
Refs. [|. 

For smaller ai, the system (||) cannot be reduced to Eq. (^, and its two-parameter family of localised solutions 
consists of three mutually coupled waves. It is interesting to note that, similar to the case of nondegenerate three- wave 
mixing Eqs. (|^) possess an exact solution. To find it, we make a substitution w = wq sech^ (ijx) , v = vasech^^ijx) 
and u — uo sech^ (77a:) , and obtain unknown parameters from the following algebraic equations 

2 9uo 2 2 

3 + 4x^0 3 

valid for rj = and a = ai = 1. Equations (|^) have two solutions corresponding to positive and negative values of 
the amplitude. This indicates a possibility of multivalued solutions, even within the class of exact solutions. 

In a general case, three- wave solitons of Eqs. (^ can be found only numerically. Figures 1(a) and 1(b) present two 
examples of solitary waves for different sets of the mismatch parameters a and ai. When ai :s> 1 [see Fig. 1(a)], 
which corresponds to an unmatched SFM process, the amplitude of the third harmonic is small, and it vanishes for 
ai ^ 00 according to the asymptotic solution of Eq. (^) discussed above. 

To summarise different types of three-wave solitary waves, in Fig. 2 we plot the dependence of the total soliton 
power defined as 

P = dx (\w\'' + 4|t.|2 -I- -\uA , (6) 

J-00 \ X J 

on the mismatch parameter ai, for fixed a — \. It is clearly seen that for some values of ai (including the exact 
solution at Q!i = 1 shown by two filled circles), there exist two different branches of three- wave solitary waves, and only 
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one of those branches approaches, for large values of ai, a family of two-wave solitons of the cascading limit (Fig. 2, 
dashed). The slope of the branches changes from negative (for small ai) to positive (for large ai), indicating a possible 
change of the soliton stability. However, the soliton stability should be defined in terms of physical parameters, and 
in the case of two-parameter solitons as we have here, the stability threshold is determined by a certain integral 
determinant condition, similar to that first derived for the three- wave mixing problem |^]. 

Ratios of the maximum amplitudes of the soliton components for the three-wave solitons of the lower branch in 
the model are presented in Fig. 3, where the upper dashed curve is the asymptotic limit of two-wave solitons for 
ai — > cxD. Soliton solutions of the second (upper) branch in Fig. 2 correspond to large values of the total power and 
they have been verified numerically to be unstable. 

The analysis of the asymptotics for Eqs. (||) suggests that localised solutions should not occur for ai < 0. However, 
we reveal the existence of an extended class of very robust localised solutions which we classify as 'quasi-solitons' pO[ , 
solitary waves with small-amplitude oscillating tails. In principle, such solitons are known in one-component models 
(see, e.g., Ref. JlT|) but here the nonvanishing tails appear only due to a resonance with the third- harmonic field [see 
Fig. 4(a)]. Such solitons are expected to be weakly unstable, and this is indeed demonstrated in Fig. 4(b) for rather 
long propagation distances. 

Existence of quasi-solitons for any value of negative phase-matching with a higher-order harmonic field indicates 
that all two-wave quadratic solitons can become unstable due to an additional SFM process. This is confirmed in 
Figs. 4(c,d) where we present the results of numerical simulations of the dynamics of an initially launched two-wave 
soliton for two cases, positive and negative phase-matching of a SFM process. For ai > [see Fig. 4(c)], a very 
small harmonic (wmax ~ 0-1) is generated and the initial two-component beam converges to a three- wave soliton. In 
contrast, for ai < [see Fig. 4(d)], the input beam decays rapidly into radiation and diffracting harmonic fields. 

In conclusion, we have investigated, analytically and numerically, multistep cascading and nonlinear beam propa- 
gation in a diffractive optical medium and introduced a novel type of three-wave parametric spatial optical solitons, 
multistep cascading solitons. The detailed analysis of the soliton stability, the effect of walk-off, higher-dimensional 
and spatio-temporal effects are possible directions of the future research. 

The authors are indebted to K. Koynov, R. Schick, and E. Kuznetsov for useful discussions. The work has been 
partially supported by the Australian Photonics Cooperative Research Centre and the Australian Research Council. 
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FIG. 1. Examples of three-wave solitary waves of the model (|^) for (a) a. = 0.05, ai — 5, and (b) a — 5, ai — 0.35. 
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FIG. 2. Two branches of multistep cascading solitons shown as the total soliton power P vs. the parameter ai for a — 1 
and X = 1- Filled circles show the analytical solutions. The lower branch approaches a family of two- wave quadratic solitons 
(for Qi — » oo) shown by a dashed line. 
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